A fast and robust approach to long-distance quantum communication with atomic 

ensembles 
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Quantum repeaters create long-distance entanglement between quantum systems while overcom- 
ing difficulties such as the attenuation of single photons in a fiber. Recently, an implementation of 
a repeater protocol based on single qubits in atomic ensembles and linear optics has been proposed 
[Nature 414, 413 (2001)]. Motivated by rapid experimental progress towards implementing that 
protocol, here we develop a more efficient scheme compatible with active purification of arbitrary 
errors. Using similar resources as the earlier protocol, our approach intrinsically purifies leakage 
out of the logical subspace and all errors within the logical subspace, leading to greatly improved 
performance in the presence of experimental inefficiencies. Our analysis indicates that our scheme 
could generate approximately one pair per 3 minutes over 1280 km distance with fidelity (F > 78%) 
sufficient to violate Bell's inequality. 



I. INTRODUCTION 

Quantum communication holds promise for the secret 
transfer of classical messages as well forming an essential 
element of quantum networks, allowing for teleportation 
of arbitrary quantum states and violations of Bell's in- 
equalities over long distances While experimental 
and even commercial implementation of simple quantum 
communication protocols are well established [2, H[ , ex- 
tending these techniques to distances much longer than 
the attenuation length of optical fiber remains a challeng- 
ing goal due to exponential attenuation of transmitted 
signals. Quantum repeaters [1,0,0] overcome the expo- 
nential time overhead associated with fiber attenuation 
and other errors by using a quantum memory and local 
quantum computation. 

Several promising avenues for quantum repeater imple- 
mentation include both atomic ensembles [7| and using 
few qubit quantum computers, such as neutral atoms in 
cavity QED [1, ion traps and solid-state single 
photon emitters (111] - Experimental progress (l2l. [l3l fl4| 
towards realization of the DLCZ protocol 0] has been es- 
pecially rapid, with many building blocks demonstrated 
in the laboratory. The experimental challenge is now 
shifting towards the realization of scalable quantum re- 
peater systems which could yield a reasonable communi- 
cation rate at continental distances (> lOOOfcm). Thus, 
the DLCZ protocol should be examined and adapted to 
practical experimental considerations, allowing to remove 
imperfections such as the finite efficiency of retrieval 
and single-photon detection and fiber length fluctua- 
tions. Our approach extends the DLCZ protocol, keeping 
the experimental simplicity of the original scheme while 
avoiding fundamental difficulties due to these expected 
experimental imperfections. 

This paper is organized as follows. In Sec. II we will 
review the DLCZ protocol and describe our new approach 
which uses a new basis to encode each qubit. Section III 
compares both the DLCZ protocol and our approach in 
the presences of imperfections. Section IV estimates the 



time scaling of our approach and compares three specific 
implementations. Section V summarizes our results. 



II. ATOMIC-ENSEMBLE-BASED QUANTUM 
REPEATERS 



A. The DLCZ protocol: a review 

The DLCZ protocol Q starts with entanglement gen- 
eration (ENG) by counting the interfering Stokes photons 
scattered from a pair of distant atomic cells x and y. This 
generates an entangled state 
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with §1 and the creation operators of spin- wave modes 
in the two cells respectively, and 4> the phase differ- 
ence between left and right channels for Stokes pho- 
tons 15[. Then, entanglement connection (ENC) is per- 
formed on two pairs of entangled atomic cells l^)^ yc 
and \£,rt> 2 ) xc y R , obtaining a further separated entangled 
pair |£0i+0 2 ) XL Vr probabilistically. The ENC step pro- 
vides built-in purification against many imperfections - 
photon loss, atomic excitation loss and dark counts. In 
the final step, post-selection is used to obtain an effec- 
tively polarization entangled state 
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from two parallel pairs 1^)^,^ \£,<t>) x ^ V2 , wmch 
comes static phase errors (time independent </>'s). 

There are two important merits of the DLCZ protocol. 
First, it has intrinsic purification of errors due to photon 
loss (in the fiber, the quantum memory, and the photon 
detector) and significantly relaxes the experimental re- 
quirement for quantum repeater. In addition, the time 
scaling of the DLCZ protocol is always sub-exponential 
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FIG. 1: Repeater components: (a) Entanglement generation 
(ENG). (b) Entanglement connection (ENC); indicated oper- 
ations: retrieve be and ac [additional 45° rotations only for 
the first level], join on polarizing beam splitter (PBS), detect 
in ± basis conditioned on one photon per output, and finally 
adjust the phase, (c) Entanglement purification (ENP); in- 
dicated operations: retrieve ai,6i and 02,62 [additional 45° 
rotations to purify phase error], interfere 01, 0,2 on PBS (same 
with 61,62), restore 03,63 conditioned on single photon at 04 
and 64 respectively, and finally adjust the phase. 



and very close to polynomial when the retrieval and de- 
tection efficiency is high. However, the DLCZ protocol 
does not purify all kinds of errors. For example, time 
dependent <p's (due to fiber length fluctuation) induce 
phase error, which cannot be taken out as a common 
factor in Eq.(j5J), since the two pairs of entangled atomic 
cells are not produced at the same time. Such phase er- 
ror is accumulated and doubled after each level of ENC. 
In addition, combined photon loss during ENG and ENC 
may also induce phase error not purified by the DLCZ 
protocol. Furthermore, the DLCZ protocol (dashed line 
in Fig. 2) still has a significant time overhead for long 
distances, because of the super-polynomial scaling in the 
presence of realistic imperfections. For instance, non- 
ideal retrieval and detection efficiency (7/ < 1) during 
ENC introduces a large vacuum component, suppresses 
the success probability of later ENC, and consequently 
slows down the protocol. 

Motivated by these issues, we will extend the DLCZ 
protocol, mitigating the above errors. 



B. New approach 

We now consider a different approach in which two 
atomic cells are used at each node a, labeled (a, H) and 
(a, V), to store one qubit, a. The qubit is defined as one 
single spin-wave excitation shared between two cells: 
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FIG. 2: Comparison between the DLCZ protocol and the 
new scheme (NS) (without active entanglement purification 
(ENP)). For each distance, we optimize over the choice of the 
control parameters (the half distance between neighboring re- 
peater stations, Lo, and the elementary pair generation prob- 
ability, p c ). With targeting fidelity F — 90%, we find the most 
efficient implementations to create the polarization entangled 
state (Eq.(f2J) for both the DLCZ (circled black dashed line) 
and the new scheme (squared blue solid line). The fiber at- 
tenuation length is L a tt — 20km, with no dynamical phase 
error. The main plot: we show the relationship between the 
(optimized) average creation time t avg and the final distance 
L for both schemes, and the empirical estimate (Eg. pip ) 
of the time scaling for the new scheme (blue dotted line). 
Over long distances (L > 320fcm), the polynomial scaling of 
the new scheme is more favorable than the super-polynomial 
scaling of the DLCZ protocol. The inset: we plot the fideli- 
ties of the intermediate distances {Lint = 160, 320, 640 and 
1280fcm), to create polarization entangled states (L = 1280, 
F — 90%), with the optimized choice of the control param- 
eters (L , Pc ) = (80,0.0027) and (40,0.0081) for the DLCZ 
(tav S ~ 1900 sec) and the new scheme (t a v g ~ 380 sec), respec- 
tively. The optimized choices of the control paramters are 
detailed in Table HVl and IV1 



When the stored spin waves are converted back into pho- 
tons, the photons have a polarization (H or V) consistent 
with that stored in the originating cell. This qubit basis 
allows projective measurements along any qubit states, 
e.g., |±) = (\H) a ± \V) a ) /\2, using linear optical oper- 
ations and photon counting [17|. We will show that in 
this logical basis it is possible to perform entanglement 
purification (ENP) @ to reduce errors within the logical 
subspace, including phase fluctuation. Since ENP can 
suppress errors within the logical subspace which occur 
with probability q to O (q 2 ), only a few ENP levels are 
needed to obtain high fidelity entanglement. 

We now describe our procedures for ENG, ENC and 
ENP. ENG (Fig. 1(a)) is similar to that of the DLCZ 
protocol, but here two parallel entangled pairs are gen- 
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crated between a and b: 

(a, if) (b.H) If*) 
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The entangled states are prepared in the quantum mem- 
ory, so no simultaneity is required for creating the two 
states comprising \ty ENG y For small excitation proba- 
bility p c , the whole generation only takes time O (l/p c ), 
in contrast to O (l/pf) for schemes requiring simultane- 
ity e.g., coupling between trapped atom and photon [Tcj 
or parametric down conversion Errors from multi- 
photon events occur only with probability pj,, and are 
considered in later analysis of imperfections. 

The first level of ENC converts two | ^ ENG ) states (one 
between cll and be, the other between ac and &r) into 
polarization entangled states \& + ) ab = \H) a \H) b + 
\V) aL \ V) bR . Only four out of the sixteen terms in 
the Schmidt decomposition of I f E)VG ) , \^ ENG ) 

r I I a L bc ' I aeon 

have any contribution to the output state; the remainder 
are eliminated by projective measurement during ENC, 
reducing the probability of success for ENC from 1/2 to 
1/8. At higher levels of ENC, the operations correspond 
to standard entanglement swapping 0, Q , where 
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leads to an entangled pair between L and R with proba- 
bility 1/2, as detailed below. 

The procedure for ENC is illustrated in Fig. 1(b). 
First, the spin waves stored in qubit be and ac are re- 
trieved into photons. At the lowest level of ENC, the 
polarization of the photons is rotated 45°. The rotations 
transforms \HV) ac{m bc) into (\HH) - \VV)) ac 
because for bosonic fields Hal 
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Thus, after the polarizing beam splitter (PBS) there will 
be at least two photons at one output. For incoming state 
I ^engx \jgENG\ a jj seven terms containing 

two excitations in at least one pair of cells in the center 
repeater node (such as \HV) ac ^ or bc ^) do not contribute 
to the click patterns with one photon at each output. 
Five terms containing two excitations in one of the left 
or right repeater nodes (e.g., \HV) a , QI bg \) have at most 
one excitation retrieved from be and ac, which is insuf- 
ficient to give two clicks. Therefore, only the four terms 
remaining can give the correct photon detector click pat- 
terns. 

For all levels of ENC, the photons are then joined on 
the middle PBS and the number of photons at two out- 
puts are counted in the {|+) , |— )} basis. With probabil- 
ity 50%, there is one photon at each output, and the con- 
nection is successful; otherwise the process is repeated. 



Operation 


Transform of $ ± ) , $ ± ) , 


Retrieve bc,a>c 
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|0000) + 1 1111} 


Transform be , ac 
at PBS 


±\0) aL (\HV) ac )m bH 
±|1>« £ (\ hv K) \°h R 


One photon per 
mode (p = 0.5) 


|0000) + |1111> 


Detect in ±, re- 
sults m,m! 


m aLbR + mm' \ll) ahbR 


Phase shift mm' 


lar.bn 



TABLE I: Entanglement connection procedure applied to 
1$ ) inputs for entangled pairs between «l and be, and ac 
and bu- For clarity, we introduce |0) = \H) and |1) = |V) 
to represent logical states (i.e. states with exactly one exci- 
tation), \HV) for non-logical states with two excitations, and 
| vac) for states with no excitation (sometimes omitted). We 
assume H photons pass through and V photons be reflected 
at the middle PBS. 



If the two photons have orthogonal polarizations, a bit 
flip a \H) +j3 \V) -»■ a \V) + \H) is applied to a L [H]. 
At higher levels of ENC, where the 45° rotations are 
not necessary, the bit flip is replaced by the phase flip 
a \H}+/3\V) a \H)-/3\V), as detailed in Tabled] [l| 

The third component is ENP (Fig. 1(c)) which obtains 
a high fidelity entangled pair from two pairs. Our proce- 
dure uses polarization entangled photons and is similar 
to recent experimental investigations . During entan- 
glement purification of bit errors (bit-ENP), the qubits 
from two parallel pairs p ai ,bi and p a2 M are retrieved from 
the quantum memory and joined at PBSs. The photons 
for two upper outputs are stored into quantum memory 
03 and 63. The photons for the lower outputs 04 and 64 
are counted in {|+) , |— )} basis. With probability 50%, 
there is exactly one photon at each lower output, and the 
purification is successful; otherwise two new pairs are cre- 
ated by restarting the process. If the two photons have 
orthogonal polarizations, a phase flip is applied to 03. An 
example of purification of bit-error is presented in Table 
ITU During purification of phase errors (phase-ENP), ad- 
ditional 45° rotations are applied to the retrieved qubits 
and the bit flip is replaced by the phase flip. The addition 
of 45° rotations effects the basis transform |$~) <-> \^ + ), 
leading to purification of errors of the other type. The 
truth table of phase-ENP is listed in Table [TTT1 Bit (or 
phase) errors can be non-linearly suppressed to the sec- 
ond order during bit-ENP (or phase-ENP) [2r| . 

The three components described above for quantum 
repeater protocol only use atomic cells, linear optics, and 
photon number counting. We remark that the duration of 
the retrieved anti-Stokes pulse can be made long (> 1/xs) 
compared to the detector recovery time by adjusting the 
intensity and duration of the retrieval pulse. This enables 
photon number counting of the anti-Stokes pulse [l9| . 
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Operation 


Transform of I*") , * n '\ 


Retrieve 
difci, a2&2 


|0001) + n |1101) + n' |0010) + nn' |1110) 


Interfere oi, a 2 on 
PBS; same with 
61,62 


]oo> aaa4 (l^> b4 ) 

+ n(\HV) a3 )\ll) b3b4 
+ (\HV) aA )\00) b3b4 
+ nn'\ll) a3ai (\HV) b3 ) 


One photon per 
lower output 






TABLE II: Bit entanglement purification acting on 
|$ + ) | , 3' + ). The same notation as Table U is used. For such 
given input states, the click pattern never matches the right 
pattern, and bit-error of single input qubit can be filtered 
completely. The remaining steps to preserve the desired Bell 
states are detection of 04, 64 in ± basis, a phase shift n ■ n' 
based on the measurement results, and storing of the travel- 
ling photons (03,63) into atomic ensembles. 
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TABLE III: Truth table for phase-ENP. Each element give 
the possible output state after the purification operation. 
("— " for cases with no outputs.) 



III. NOISE AND IMPERFECTIONS 

A. Non-logical errors 

We now examine the performance of our new scheme 
by considering the role of errors, starting with how im- 
perfections due to inefficiency limit the protocols. Pri- 
marily, we find that inefficiency takes logical states into 
two types of non-logical states - those with too few exci- 
tations (vacuum type) and those with too many excita- 
tions (multi- excitation type). We represent these errors 
by density matrix ir vac (a mixed state with at most one 
excitation between both pairs of cells) or Tt m uiti (a mixed 
state with at least one pair of cells with more than one 
excitation). The normalized density matrix after mth 
ENC (Fig. 1(b)) can be written as 

(m) (m) (m) ( m ) ( m ) (to) (to) 

^a L ,b R ~ Plogicr logic T Pvac "vac ~ P multi" multi \° > 

where the m-dependent operator p[™] c is the density ma- 
trix within logical subspace; Kiac and T^nuiu also depend 

on m; and the coefficients f>j™] c , Pi2l and p^liu are tne 
probabilities for the logical, vacuum and multi-excitation 
types, respectively. 

After the first level of ENC, piW ~ 1 



p c < 1. We can demonstrate that these three prob- 
abilities remain stable (see Appendix B) for all higher 
levels of ENC, by considering the un-normalized state 
after (to + l)th ENC for the new approach: = 
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' TpWlogicP multi V 1 


+ O (Vc)) 



where the probability for the new logical error from the 
multi-excitation states (accompanied by photon loss) is 



V, 



(m+l 



L ~ (1 - V)Pc 



The total logical error probability pi™ has two contri- 
butions: the accumulated logical errors from both input 



pairs for ENC, pir), and the new logical error 

p£+V-2p£)~(l-r,)p c . 
We can calculate for the new scheme, 



(6) 



(7) 



V and p^ ulti 



A more detailed calculation (see Appendix B) , in 
which ir vac and n mu iu are further divided into sub- 
spaces with different number of excitations (e.g. Tr vac is 
subdivided into zero-excitation and one-excitation sub- 
spaces), verifies the stability of the probability distribu- 
tion of P;™] c , Pv"o!c and p^2 la - (Similarly, dark count can 
also induce errors in logical subspace with probability 
~ Pdark (1 ~ i) s ), which is however negligible due to very 
low dark count probability p c iark) 

For the DLCZ protocol, only two cells are used to store 
entanglement. Besides the logical states (single excita- 
tion in two cells), we can similarly define the vacuum 
states (with no excitation) and multi-excitation states 
(with two or more excitations). Contrary to our ap- 
proach, the probability distribution is not stable - both 
vacuum and multi-excitation probabilities increases with 
distance (see Appendix A). The vacuum probability soon 
becomes the dominant term, which reduces the suc- 
cess probability of ENC significantly, resulting in super- 
polynomial (but still sub-exponential) time scaling (Fig. 
2(b)). The logical error probability for the DLCZ pro- 
tocol has the same form as Eq.© up to the coefficient, 
but the ratio P ™u\ulP logic (thus pi™J) grows with distance 
(see Appendix A), which accounts for the sharp decrease 
of fidelity for the DLCZ protocol in Fig. 2(a). To main- 
tain good final fidelity, the initial error p^ uHi (and p c ) 
should be very small, which demands much longer gener- 
ation time of an elementary pair for the DLCZ protocol 
(Fig. 2(b)). 

In essence, by requiring at least one excitation in the 
ensemble, our qubit subspace is automatically purified of 
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vacuum and multi-excitation type errors during ENC. 
The closest analog to the DLCZ protocol is our new 
scheme without ENP, i.e., only ENC. At longer distances, 
our approach is further improved in comparison to the 
DLCZ protocol due to the reduced amplitude of vacuum 
terms. 



B. Logical Errors 

So far, we have only considered the effects of ineffi- 
ciency that maps states between logical and non-logical 
subspace by changing the number of excitations. Besides 
inefficiency, there are other imperfections, which pre- 
serve the number of excitations but induce errors within 
the logical subspace, such as interferometric pathlcngth 
fluctuation and linear optical misalignment. For exam- 
ple, the interferometric pathlength fluctuation leads to a 
phase difference 5, which changes \ 1 & ENG ) a b into a mix- 
ture: 

\* ENG ) a , b - ^ (\ H )a \V) b + e* \V) a \H) b ) (8) 
+ | J ff^)Jvac) h + e 2 ^|vac) a |W) h . 

where (j) is static phase difference between left and right 
channels. Since the last two terms with \HV) will be 
removed during the first level of ENC, the static phase 
4> has no effect. However, the probability of being in 
an undesired logical state Y$>~) is sin 2 |. The first 
level of ENC with the combined two inputs of "f - and 
ty + gives $ _ , producing a phase error with probability 
%Pphase—err — 2 (sin 2 |), proportional to the variance of 
the interferometric phase fluctuation. This error will be 
amplified during subsequent ENC's, because the survival 
probability of the state <I> - (the logical error) is twice as 
much as that of <J> + (the desired component). 

We expect that there is little correlation in phase fluc- 
tuation between different sections of the fiber, and the 
variance of the phase fluctuation is proportional to the 
length of the fiber 

(S 2 ) = 2DL 0l (9) 

where D is the phase diffusion coefficient of the fiber. If 
the phase fluctuation satisfies gaussian distribution, the 
phase error probability 




For example, D = 10 radr /km, Lq = 10km, and 

Pphase—err ~ 0.5%. 

Also, a small probability (p err ) of linear optical mis- 
alignment per ENC or ENP step is modeled as depo- 
larizing errors. Later, we will demonstrate that errors 
within the logical subspace restrict the final fidelity of 
the DLCZ protocol, while for our new approach addi- 
tional active purification can correct such logical errors 
to achieve high fidelity. 



IV. SCALING AND TIME OVERHEAD FOR 
QUANTUM REPEATER 

A. Scaling analysis 

Based on the calculation of the success probability at 
each level of connection/purification, we can obtain the 
estimated average time for various schemes. In Fig. 2, 
we compare our approach to the DLCZ protocol. For the 
DLCZ protocol, the average creation time for a distant 
pair contains a super-polynomial contribution (but still 
sub-exponentially) with distance, due to instability of the 
vacuum component. For our new scheme, the scaling is 
strictly polynomial with distance, t avg oc L Q , where the 
exponent a — a (77) explicitly depends on the efficiency. 

We remark that the DLCZ protocol is slightly more ef- 
ficient for short final distances (L < 160fcm). The DLCZ 
protocol can skip the entanglement connection and ex- 
ploits post-selection to create the polarization entangled 
state, while the new scheme requires the first level of 
ENC to eliminate those unwanted components and pre- 
pare the polarization entangled state. The post-selection 
has success probability 1/2, two times more efficient than 
the success probability 1/4 for the first level of ENC of 
the new scheme. 

For the new scheme (without ENP), we can use the 
stable probability distribution to estimate the average 
time: 

t^^toiL/Lof^i 1 - 5 ^^) (11) 

where to = ]Trj^ c a ' eLo ^ La " i s the elementary pair genera- 
tion time, pc is the elementary pair generation probabil- 
ity, Lq is half the distance between neighboring repeater 
stations, and L is the final distance. The exponent can be 
understood as an overhead from the finite success prob- 
ability for ENC, penc ~ n l,^Z 2 $ ■ Thc constant 1.5 in 

z{A rj) 

Eq. pT|) is the empirical estimate of the overhead from 
the waiting time to obtain two independent pairs ver- 
sus the single pair. In Fig. 2, the differences between 
the simulated data and Eq. pT|) are attributed to the em- 
pirical factor (1.5) and an overall factor from the over- 
estimate of the success probability for the first level of 
ENC. According to Eq.([7]), one can always reach good 
final fidelity if the elementary pair generation probabil- 
ity p c scales as Lq/ L. Therefore, the average distant pair 
generation time scales exactly polynomially with distance 

t avg oc L a , with a = 1 + log 2 (1.5) + log 2 ( |^^) ) ■ 

B. Comparison between different schemes 

Besides the DLCZ protocol and our new approach 
without ENP, we now consider a scheme with ENP, which 
has one phase-ENP after the second level of ENC. We 
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FIG. 3: Average time versus final fidelity with no phase er- 
rors. For fixed final distance L — 1280, we plot the (final) 
fidelity dependence of the (optimized) average pair creation 
time, for the DLCZ protocol (black dashed lines), new scheme 
without ENP (blue solid lines), and new scheme with ENP 
(red dashdotted lines), with efficiency rj to be 90% (circles), 
and 95% (no circles). 



FIG. 4: Average time versus final fidelity with phase errors. 
This shows the same optimized approaches as Fig. 3 but with 
finite interferometric pathlength fluctuation characterized by 
the phase diffusion coefficient D = 10" 3 r ad 2 /km. Inclusion 
of ENP (red dashdotted lines) yields a dramatic improvement 
in time for high fidelity operations. 



may compare these schemes by using t — F plots - a 
parametric plot of t avg and F as a function of excitation 
probability p c . For given noise model and efficiency 77, a 
repeater scheme corresponds to a curve on t — F plane. 

In the absence of interferometric pathlength fluctua- 
tion (Fig. 3), the new approach without ENP is about 
5 times faster than the DLCZ protocol, for ?y = 90%. 
As given by the previous discussion, this improvement 
is due to better control of inefficiency-induced imperfec- 
tions. There is a time overhead for the new approach 
(with ENP) as compared to the new approach (without 
ENP). Within each implementation, the higher the effi- 
ciency 77, the faster the quantum repeater. For high final 
fidelity (1 — F < 10%), the curves approach straight lines 
with slope —1, because t oc p// 1 oc (1 — F) 

When interferometric pathlength fluctuation (leading 
to initial phase error) is non-negligible, active ENP is 
needed. We use a diffusion model for the pathlength fluc- 
tuation, as detailed in Sec. IIII Bl In Fig. 4, t — F curves 
are plotted, assuming the phase diffusion coefficient D = 
I0~ 3 rad 2 /km, corresponding to p p hase-err ~ 0.5% over 
Lq = lOfcm. Unlike Fig. 3 where only inefficiency is con- 
sidered, there is an upper bound in final fidelity for each 
implementation. Both the DLCZ protocol and new ap- 
proach (without ENP) suffer from the initial phase error, 
with final fidelity no more than 65%, while new approach 
with ENP maintains high final fidelity up to 97%. For 
high retrieval and detection efficiency (?? = 95%), new ap- 
proach (with ENP) can produce 1280fcm entangled pairs 
with fidelity 90% at a rate of 2 pairs/hour, even in the 
presence of substantial dynamical phase errors. 



V. OUTLOOK 

In summary, our new approach to long distance quan- 
tum communication uses a different qubit basis which 
prevents the growth of vacuum and multi-excitation 
probabilities. This keeps the ENC success probability 
high and error probability low, and leads to true polyno- 
mial scaling even in the presence of realistic inefficiencies. 
We can achieve a bandwidth of 1 (or 2) entangled pair(s) 
per six minutes for 77 ss 90% (or 78%) and negligible ini- 
tial phase error. The new approach also allows active 
entanglement purification, which combined with built-in 
purification of transmission loss errors allows purification 
of arbitrary errors in quantum communication. 

Although the present approach shows a dramatic im- 
provement in communication rates and robustness com- 
pared to original the DLCZ protocol, the bandwidth 
remains relatively slow, even when very high efficien- 
cies and very long-lived quantum memory are assumed. 
While such hi gh e fficiencies might ultimately be achiev- 
able (see Ref. [2l| for recent progress), other approaches 
need to be considered that can further improve the effec- 
tive communication bandwidth. For example, we can use 
many cells per node to improve the bandwidth. In this 
case, the improvement is at least linear with the num- 
ber of cells, making it possible to realize long distance 
(1280fcm) entangled state generation bandwidth of the 
order of one pair per second. 

A simple Monte Carlo optimization of efficient use of 
the cells shows that we can increase the bandwidth by a 
factor of r = M 112 , where M is the increment of factor of 
physical resources. Recently, it has also been suggested 
that multiple cells can be used to further facilliatc the 
quantum repeater, in the presence of memory errors from 
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L (km) 


20 


40 


80 


160 


320 


640 


1280 


2560 


5120 


10240 


tavg (sec) 


0.0013 


0.0065 


0.046 


0.44 


5.8 


95 


1900 


5.7 x 10 4 


3.0 x 10 6 


3.0 x 10 s 


Lo (km) 


10 


20 


40 


40 


40 


80 


80 


80 


80 


80 


Pc 


0.26 


0.18 


0.13 


0.040 


0.011 


0.010 


2.7 x 10~ 3 


7.1 x 10~ 4 


1.8 x 10~ 4 


4.6 x 10~ 5 


Ffin 


90% 


90% 


90% 


90% 


90% 


90% 


90% 


90% 


90% 


90% 



TABLE IV: Minimized average time tavg and optimized control parameters (Lo , p c ) for the DLCZ protocol. This table provides 
detail information for Fig. 2. 



L (km) 


20 


40 


80 


160 


320 


640 


1280 2560 5120 10240 


tavg (sec) 


0.0051 


0.020 


0.10 


0.68 


5.4 


45 


380 3.3 x 10 3 2.9 x 10 4 2.6 x 10 5 


Lq (km) 


5 


10 


20 


40 


40 


40 


40 40 40 40 


Pc 


0.26 


0.17 


0.11 


0.087 


0.037 


0.017 


8.1 x 10" 3 4.0 x 10" 3 2.0 x 10~ 3 9.7 x 10~ 4 


Ff in 


90% 


90% 


90% 


90% 


90% 


90% 


90% 90% 90% 90% 


TABLE V 


: Miminized 


average 


time tavg 


and optimi2 


ed control parameters 


(Lo , p c ) for the new scheme. This table provides 



detail information for Fig. 2. 



the quantum memory [22l ]. 
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APPENDIX 

We now justify two main claims used in the previous 
discussion. (1) For the DLCZ protocol, the probabil- 
ity ratios pi™l/p ( ™J ic and P { ,nuiti/p\7 g lc increase with the 
nesting level m. (2) For the new scheme (without ENP), 

the probability ratios piai/^Sc and vZliulvZgXc remain 
almost independent of the nesting level m. The states 
of logical, vacuum and multi-excitation types depend on 
the repeater protocol. For example, the entangled logical 
state for the DLCZ protocol has total of one excitation 
stored in two remotely entangled cells (e.g. Eq.((T])), while 
the entangled logical state for the new scheme has two 
excitations stored in four cells, with one excitation and 
two cells on each side (e.g. Eq.@). Thus, the defini- 
tions of vio g \^ Pvac and p^ Hi are different for the two 
schemes. 



APPENDIX A: NON-LOGICAL STATES FOR 
THE DLCZ PROTOCOL 

We start with the DLCZ protocol. First, we decom- 
pose the density matrix (for a pair of distant atomic cells 
x and y) into components with different excitation pat- 
terns, neglecting the inter-pattern coherence 

Px,y = POOTOO + PlOTTlO + Pll^l! + P207T20 H (Al) 



Note added. - Since completion of this work, a preprint 
by Chen et al. (23| describing a similar approach to im- 
prove the DLCZ protocol has appeared. 



I 

where Pij^ij is the projected density matrix to the 
subspacc spanned by the Fock states of the cells 

\j)y i \f)x l*)y}' W1 th probability py and normalized 
density matrix itij. For the DLCZ protocol, the vacuum 
type of state is 7Too, the logical type of state is ttiq, and 
the rest belong to the multi-excitation type. We may also 
introduce the notation corresponding to Eq.([5]) 

Px,y = PlogicPlogic + Pvac^vac + Pmulti^multi 

(A2) 

withp„ ac = poo, piogic = Pio, and p mu iti = Pn+P2oH • 

Since we are only interested in the coherence proper- 
ties for the logical type 7Tio, we only keep track of the 
probabilities for the vacuum and multi-excitation types 
and neglect their coherences. From the symmetry for the 
two cells, we have 

Ttf = 71* = l ~ (\i) x (i\ ® \j) y (j| + \j) x (j\ (» \i) y . 

(A3) 

In the rest of this subsection, we will add a superscript to 
■Kij only when we want to keep track of the coherence for 

that specific term. For example, 7r^^ indicates that there 
is coherence between the states 1 1)^ |0) and (0)^ |1) , 
after the mth level of ENC. 

Since the DLCZ protocol requires that the probabil- 
ity for the multi-excitation states should always be much 
smaller than the probability for the logical states (other- 
wise, a large fraction of multi-excitation states, accompa- 
nied by photon loss, can induce significant logical errors), 

, (m) (m) I (™)\ ; i i • 

we regard p n , p 2a (<C p\ ) as perturbations. 
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We denote the entangled state after the mth ENC as 



,(m) 



(m) 

Poo 7roo+Pio w w +Pn nn+Pm ^20+- • ■ (A4) 



(m) (m) (m) 



(m) 



and 



We now connect two such entangled states p 
Pxr]vr f° r the (m + l)th ENC via the superoperator 

&ENC 



~(m+l) 
rx L ,y R 



£enc 



(A5) 



where pi™^ is the unnormalized density matrix for the 
entangled state after the (m + l)th ENC. Since £enc is 
a linear operator and two inputs have the same state, 



£enc 



/3 



(A6) 

where Wa^p] = \£enc [tTq,^] + \£enc ka, n a ]. 

ENC [nenp] for 71^,7173 G 



'.EJVC 

Now we calculate £ sym 



(m) 

7r oo,7i'io jTriifTaO)' 



|. For example, 



psym 
°ENC 



c'Syvri 
C ENC 



psyra 
C ENC 



(m) (m) 
^10 ! w 10 

(m) 

TTlO '"""00 

(to) 
^lO > ""11 



7? (m+1)/ 7/(1 -r/) 
= 2^10 + 2 71-00 ' ^ 



fsym 
°ENC 



(m) 

^lO '^20 



7T00 

7Tn +77(1-77) 7Tio 
77 (1 - 7?) 



77 

7^0 



, 377 (1 - 77)^ 



-7T1Q 



csyva r i 

L ENC F00i ^OOj — U 
csyra r 1 

^_e/vc Foo, Tn J = 777rio 

^BJVC fa, 7T20] = V (1 - 7?) 7T0Q. 



(A7b) 
(A7c) 
(A7d) 

(A7e) 

(A7f) 
(A7g) 
(A7h) 



The logical type of state after the (m + l)th ENC, 

7rj™ +1 \ is the average of the states 7rj™ ' and 7Tio, 

with relative weights fp^Pio^ and 2?7 (1 _ '?) Pio^Pn + 

77 (1 - rj)PwP20 + 2 VPoo l) Pii \ respectively. 

Then we calculate the (unnormalized) density matrix 
after the (m + l)th ENC, 



to the leading order with respect to the perturbations of 

(m+l) 



and ;4™ } 



Finally, we divide these probabilities by p\ , and 
obtain 



~(m+l) 

Poo 


/ (m 
~ 9 (^00 

" (m 


~(m+l) 

P10 


VPio 


-(m+l) 
Pil 


(m) 
~ (m) 

Pio 


-(m+l) 

P10 


~(m+l) 
P20 


(m) 
^ P20 


-(m+l) 

P10 


~ M ■ 

P10 



(1-'/) 



(m) 

Poo 

(m) 
PlO 



(AlOa) 
(AlOb) 
(AlOc) 



Since the normalization does not change the relative 
ratio between the probabilities, the above purterbative 
estimate tells us that the fractions for both the vacuum 
state (7100) and the multi-excitation state (%n) are at 

least doubled, relative to the logical state 7rJ™ +1 '', 
each ENC. 

In terms of the notation corresponding to 

(Pvac = P00, Plogic = PlO, and Pmulti = Pll + P20 

we have 



after 

Eq.© 

f ■■•), 



Pvac 

(m+l) 
Plogic 

~(m+l) 
Pmulti 

-(m+l) 
Plogic 



-(m+l) 

Poo 

-(m+l) 
PlO 

(m) 
O Pmulti 

(m) 
Plogic 



> 2 



(m) 
Pvac 



Pi 



ogic 



(Alia) 
(Allb) 



The ratio of non-logical states to logical states is at least 
doubling with distance. As discussed in the main text, 
it is these unstable non-logical states that leads to the 
super-polynomial scaling for the DLCZ protocol. 



APPENDIX B: NON-LOGICAL STATES FOR 
THE NEW SCHEME 

For the new scheme, we can similiarly decompose the 
density matrix (following Eq.©) 

P = POOTTOO + PlOTTlO + PllTTll + P20T20 

+ P217T21 +P30130 ' ' ' (B12) 



~ (m+l) _ 5 (" 1 + 1 ) 7rnn + 

"xL,yn — POO 00 T Pio n 10 

-(m+l) 



Pll 



, -(m+l) 
7T11+P20 ^20 



with the probability coefficients 

-(m+l) (m) / (m) (1 — 7?) (m) 

Poo ~ Wio Poo + — o — PlO 



~(m+l) 7) ( m ) ( m ) 

Pio ~ 2 p i° Pl ° 

-(m+l) ~ (m) (m+l) 

Pll ~ WlO Pll 

~(m+l) 7y (to) (m+l) 

P20 ~ 2 Pl ° P20 



(A8) 

(A9a) 

(A9b) 
(A9c) 
(A9d) 



where PijKij is the projected density matrix to the sub- 
space spanned the Fock states with i (or j) photons in a 
cell-pair and j (or i) photons in b cell-pair, with probabil- 
ity Pij and normalized density matrix 7Ty . For the new 
scheme, the vacuum type of states consists of 7Too and 
7Tio, the logical type of state is 7Tn, and the rest belong 
to the multi-excitation type. 

We use a perturbative approach, by assuming poo^ooj 
Piottio, and pniru are the dominant terms, and the rest 
terms are perturbations of order p c (terms not listed are 
of order p 2 ). We eliminate those irrelevant perturbation 
terms (e.g. j?3o7T3o), because after one level of ENC, they 
are suppressed to O (pf) . 
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Suppose the entangled state after the mth ENC is 



Pa,b =P00 n 00 + P W 7T10 
. fml 

""20, 



(m 
Pll 



i ( m ) 
+ P 2 0,|| 7r 20,|| 



(m 

P20, 



(m) 
^21,11^21,11 



(m) 
■P21.J 



T21, 



(B13) 



Notice that we need to distinguish two possible types 
of states for 7120 because they behave differently during 
ENC. The first type of states (denoted as 7r 2 o, || ) has both 
photons stored in the same cell, and after retrieval the 
photons will have the same polarization, follow the same 
path way, and trigger the photon detector(s) on the same 
side of the PBS (Fig. 2(b)). The second type of states 
(denoted as 7r 2 o,±) has two photons stored in different 
cells, and after the retrieval the photons will have orthog- 
onal polarization, split at the PBS, and trigger photon 
detectors on both sides of the PBS (Fig. 2(b)). Thus, the 
second type of states are more likely to give the correct 
click pattern and thus propagate the error to the next 
level of ENC. Similarly, we introduce ^2i,\\ an d 7T2i,_i_. 

In the rest of the discussion, we still follow the con- 
vention from the previous subsection that 71^- = nji and 
we will add a superscript m to 7r,y only when we want to 
keep track of the coherence for that specific term. 

We now connect two such entangled states p„ j, and 



p ( Z\b R for thc ( m + l ) th ENC 



~(m+l) c 
Pa L ,b R -ZENC 



(m) (m) 
Pa L ,bc ' Pac,ba 



(B14) 



where p ~t is the unnormalized density matrix for the 
entangled state after the (m + l)th ENC. 



Now we calculate ^£jvc Ki^l 



2&ENC [""as, T/3] + \&ENC Kg,7T a ] for 



» 



|too,tio 
ample, 



ir a , irp G 
For ex- 



r>sym r 
E NC POO, 7T21 



'ENC F00,7I"21,iJ = Y^IO 



(B15j) 
(B15k) 



The logical type of state after the (m+l)th 
ENC, 7rj™ +1 ' ) , is thc average of the states 7rJ™ +1 ^ 
and 7Tn, with relative weights \Pi\ P11 and 
V 2 (1 - v)Pn ] Es=||,xj4m + fPio'VSi: respectively. 

Then we calculate the (unnormalized) density matrix 
after the (m + l)th ENC, 



-(m+1) 
Pa L ,b R 



~(m+l) 



-(to+1) 



Poo ^00 + Pw 7T10 + Pi 



~(m+l) (m+1) 



'11 



. ~(m+l) . (m+1) . (m+1) 

+ ^20,11 ^20,11 + P 2 0,_L ; 7I"20,±+P 21) || '7r2l,|| 



+ P 2 1.+ ^21,+ 



with the probability coefficients 



(B16) 



~(m+l) , 
Poo ; 


^ V 2 (m) (m) 

a Y p i° Pl ° 


(B17a) 


~(m+l) , 

P10 f 




(B17b) 


-(m+1) , 
Pll ; 


r, 2 / pfrW^+J&'l&J. 

~ 2 i v +2(i- ?? )^T ) (^j,+^:!i 


,) 

(B17c) 


~(m+l) , 
P 2 0,<5 


^ ?7 2 (m) (m) 
~ ^ PlO P21,8 


(B17d) 


-(m+1) , 
Pil,6 f 


_ V 2 (m) (m) 
~ 2 ^ n ^21, 5 


(B17e) 



Finally, we divide these probabilities by p^ +1 \ and to 
the leading order (i.e. the zeroth order of p c ) we obtain 



psym 
C ENC 



(m) (m) 
"11 ' ^ll 



psym 
C ENC 



r>sym 
Z'ENC 

C ENC 



(m) 
TTii ',7T10 

(m) 
7Tu ,7T o 

(m) 
7Tu ,7r 20j( 5 



psym 
C ENC 



(m) 
Til , T21/ 



fp'Tr. [tTiq, 7Tio] 



■■ENC 
?sym 



rsym r 1 
£ _EJVC [TlO, T00J 

csym r 1 
^iJJVC LT10,7T21,||J 

r>sym r i 
£ _EATC F10) T21.J.J 



V (m + 1)/ 



= 



rf (1 - 77)7110 
,2 





v 2 (1 - v) 



v 2 

T 771,1 



TlO + ¥ 

^ 2 (1 - »?) 



T20, 



"TlO 



(B15a) 

(B15b) 
(B15c) 
(B15d) 



fT 2M + !^l 7 rii (B15e) 
r? 2 

ir^oo 



(B15f) 
(B15g) 

(B15h) 
ry 2 

yTao.i 
(B15i) 



-(m+1) 

Poo 




-(m+1) 
Pll 




-(m + 1) 

P10 




-(m+1) 
Pll 


" (m) 
Pll 


-(m+1) 
P20.S 


1 (m) (m) 
^ IPlO ^21,5 


-(m+1) 
Pll 


O (m) (m) 
2 Pll Pll 


-(m+1) 
P2M 


(m) 
. P21.5 


-(m+1) 
Pll 


" (m) 
Pll 



(B18) 
(B19) 
(B20) 
(B21) 



Furthermore, all these ratios remain constant (to order 
p c ), which justifies the claim that the probabilities for 
different types of states remain stable for all higher levels 
of ENC. 



If we further introduce plac +1 ^ 



-(m+1) _ -(m+1). 



Plosic ' =Pll and Pmulti 



^(m+1) 



-(m+1) ~(m+l) 
Poo "r PlO ' 

V d m+1) + n {m+1) 

l^6=\\,±P20,S + P 2 1,5 ' 
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we have 



& « { 1 + P&JS. + O (Be)) (B22a) 



fopac 2 Logic^ iog% 

.,2 



Ptc +1) « ^p!2cPSS(1 + foe)) 
2 

~(m+l) ?7_ (m) (m) Q , ,s 

where pi^^L, is the probability for the new logical error 



(B22b) 
(B22c) 



from the multi-excitation states (accompanied by photon 
loss) 



V 



(m+l) 
err,new 



(m) (m) 
Pil 7>11 



(1 - 7?)p c . 



(B23) 
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